An efficient algorithm to recognize local Clifford equivalence of graph states 
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In [Phys. Rev. A 69, 022316 (2004)] we presented a description of the action of local Clifford 
operations on graph states in terms of a graph transformation rule, known in graph theory as local 
complementation. It was shown that two graph states are equivalent under the local Clifford group if 
and only if there exists a sequence of local complementations which relates their associated graphs. In 
this short note we report the existence of a polynomial time algorithm, published in [Combinatorica 
11 (4), 315 (1991)], which decides whether two given graphs are related by a sequence of local 
complementations. Hence an efficient algorithm to detect local Clifford equivalence of graph states 
is obtained. 

PACS numbers: 03.67.-a 



Graph states have been studied extensively and have 
been employed in a number of applications in quantum 
information theory and quantum computing (see e.g. 
IH HI El 01 This is mainly due to the fact that these 
states can be described in a relatively transparent way, 
while they maintain a sufficiently rich structure. In this 
note we consider the problem of recognizing local Clif- 
ford (LC) equivalence between graph states. This issue is 
of natural importance in multipartite entanglement the- 
ory 0, [j| and in the development of the one-way quan- 
tum computer, which is a universal measurement-based 
model of quantum computation In the following, we 
present an efficient algorithm which recognizes whether 
two given graph states are LC-equivalent. At the heart 
of this algorithm lies an earlier result [B| of ours, which 
is a translation of the action of local Clifford operations 
on graph states in terms of a graph transformation rule 
known in graph theory as local complementation. It was 
shown that two graph states are equivalent under the lo- 
cal Clifford group if and only if there exists a sequence 
of local complementations which relates their associated 
graphs. As it turns out, local graph complementation is 
well known in graph theory (see e.g. |(| and references 
within). What is more, in ref. Q a polynomial time algo- 
rithm is derived which detects whether two given graphs 
are related by a sequence of local complementations. This 
yields for our purposes an efficient algorithm which recog- 
nizes LC-equivalence of graph states. We repeat this al- 
gorithm below. Note that the present result immediately 
yields an efficient algorithm to recognize LC-equivalence 
of all stabilizer states (and not just the subclass of graph 
states). Indeed, it is well known that any stabilizer states 
is LC-equivalent to a graph state Moreover, if a 

particular stabilizer state is given then an LC-equivalent 
graph state can be found in polynomial time, as the typ- 
ical existing algorithms used to produce this graph state 
essentially use pivoting methods, which can be imple- 
mented efficiently. 
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Before presenting the algorithm, we state some defi- 
nitions and introduce some notations. Graph states are 
special cases of stabilizer states. They are defined as 
follows: let G be a simple graph on n vertices with adja- 
cency matrix 9 [Til ] and define n commuting correlation 
operators Kj (j — 1, . . . , n), which act on Cf ™, by 

n „ 

iv 4 j II (^) *. 

fe=l 

Here <Jx , Oy , ai^ are the Pauli matrices which act on 
the ith copy of C2. The graph state \G) is the unique 
eigenvector (up to an overall phase) with eigenvalue one 
of the n operators Kj. The Clifford group C± on one 
qubit is the group of all 2 x 2 unitary operators which 
map a u to auc^iu) under conjugation, where u — x,y, z, 
for some a u = ±1 and some permutation n of {x,y,z}. 
The n-qubit local Clifford group C l n is the n-fold tensor 
product of Ci with itself. The action of local Clifford 
operations on graph states can be translated elegantly 
in terms of graph transformations. In ref. [f| we give 
such a translation and show that two states \G), \G') are 
LC-equivalent if and only if G' can be obtained from G 
by a finite sequence of local complementations. The local 
complement gi(G) of a graph G at one of its vertices 
i G V is defined by its adjacency matrix gi{9) as follows: 

g l (9) = 9 + 9 l 9f + A, (1) 

where 9 is the adjacency matrix of G, 9i is its ith column 
and A is a diagonal matrix such as to yield zeros on the 
diagonal of gi(9). Addition in is to be performed 
modulo two. In graph theoretical terminology, gi(G) is 
obtained by replacing the subgraph of G induced by the 
neighborhood of i by its complement. 

It is well known (see e.g. [l0|]) that the stabilizer for- 
malism has an equivalent formulation in terms of binary 
linear algebra. In this binary formulation, a stabilizer 
state on n-qubits corresponds to an n-dimensional self- 
dual linear subspace of Fj™. Here F2 is the finite field of 
two elements (0 and 1), where arithmetics are performed 



modulo 2. The self-duality of the subspace is with re- 
spect to a symplectic inner product < •, • > defined by 



Pv, where u, v £ F:f l and 



P = 



/ 
/ 



The binary stabilizer subspace is usually presented in 
terms of a full rank 2n x n generator matrix S, the 
columns of which form a basis of the subspace. This gen- 
erator matrix satisfies S T PS = from the self-duality of 
the space. The entire binary stabilizer space, which we 
denote by Cs, is the column space of S. It can easily be 
shown that a graph state with adjacency matrix 6 has a 
generator matrix 



S 



In the binary framework, local Clifford operations U £ C n 
correspond to nonsingular In x 2n binary matrices Q of 
the block form 



Q = 



A B 
C D 



where the nx n blocks A,B,C,D are diagonal |]|. We 
denote the diagonal entries of A, B, C, D, respectively, by 
di, bi, c,-, di, respectively. The n submatrices 



Q 



(i) — 



b, 
d, 



correspond to the tensor factors of U. It follows that 
each of the matrices Q^ 1 ' is invertible. Equivalently, the 
determinants a^di + biCi are equal to one. We denote the 
group of all such Q by C l n . 

We are now in a position to state the algorithm. Let 
\G), \G') be two states with adjacency matrices 6, 8', 
respectively, and generator matrices 



S := 



, S> := 



respectively. Then \G) and \G') are LC-equivalent if and 
only if there exists Q £ C l n such that Cqs = Cs> ■ Equiv- 
alently, this occurs iff there exists an invertible n x n 
matrix R over F2 such that 



QSR = S'. 



(2) 



If 9 and 8' are given, is a matrix equation in the 
unknowns Q and R. Note that we can get rid of the 
unknown R, as (J5J) is equivalent to 



S T Q T PS' = 0. 



(3) 



Indeed, expresses that u 1 Pv = for every u £ Cqs 
and v £ Cs>, which implies that Cqs and Cs> are each 
other's symplectic orthogonal complement. These spaces 



must therefore be equal, as any n-dimensional binary sta- 
bilizer space is its own symplectic dual, and is ob- 
tained. More explicitly, is the system of n 2 linear 
equations 

(j2 Mi*c^ + e ika k + 0' jk d 3 + Sjkbj = 0, (4) 

for all j, k = 1, . . . , n, where the 4n unknowns <Zj, bi, Cj, di 
must satisfy the quadratic constraints 



aidi + b t Ci = 1. 



(5) 



The set V of solutions to the linear equations , with 
disregard of the constraints, is a linear subspace of Ff™. A 
basis B = {61, . . . , bd} of V can be calculated efficiently 
in 0(n 4 ) time by standard Gauss elimination over F2. 
Then we can search the space V for a vector which sat- 
isfies the constraints |J5J. As @ is for large n a highly 
overdetermined system of equations, the space V is typi- 
cally low-dimensional. Therefore, in the majority of cases 
this method gives a quick response. Nevertheless, in gen- 
eral one cannot exclude that the dimension of V is of 
order 0{n) and therefore the overall complexity of this 
approach is nonpolynomial. However, it was shown in 
that it is sufficient to enumerate a specified subset V'CV 
with |V'| = 0(n 2 ) in order to find a solution which satis- 
fies the constraints, if such a solution exists. Indeed, the 
following lemma holds: 

Lemma 1 7] If dim(V) > 4, then the system Op-JSJ) 
of linear equations plus constraints has a solution if and 
only if the set 

V := {b + b' I 6,6' £ B} C V 

contains a vector which satisfies the constraints. 

The proof of lemma 1 is involved and makes extensive 
use of local graph complementation. The reader is re- 
ferred to ref. |7j for more details. Lemma 1 shows that, if 
a solution to Q-© exists, this solution can be found by 
enumerating either all |V| < 16 elements of V if dim(V) < 
4 or the 0(n 2 ) elements of V' if dim(V) > 4 and checking 
these vectors against the constraints J5J. Hence, a poly- 
nomial time algorithm to check the solvability of (@J-© 
is obtained. The overall complexity of the algorithm is 
C(n 4 ). Note that, whenever LC-equivalence occurs, this 
algorithm provides an explicit local unitary operator in 
the Clifford group which maps the one state to the other, 
as a solution (cii, 61, Ci, c?i, . . . , a ni 6 n , c n , d n ) to ©-© 
immediately yields an operator Q £ C l . 

In conclusion, we have presented an algorithm of 
polynomial complexity which detects whether two given 
graph states are equivalent under the local Clifford group. 
This algorithm leans heavily on a former result of ours, 
which is a description of the action of local Clifford op- 
erations on graph states in terms of local graph comple- 
mentation. Whenever equivalence of two graph states is 



recognized, the algorithm provides an explicit local Clif- 
ford operator which maps the one state to the other. 
Moreover, together with existing algorithms, this result 
yields an efficient algorithm which recognizes local Clif- 
ford equivalence of all stabilizer states. 
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